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DESCRIPTION 

VISUALIZATION PROCESSING SYSTEM, VISUALIZATION PROCESSING METHOD, 
AND VISUALIZATION PROCESSING PROGRAM 

5 

Field of Art 

The present invention relates to a visualization processing system, a visualization 
processing method, and a visualization processing program, and particularly, to a visualization 
processing system, a visualization processing method, and a visualization processing program 

10 which are adapted for visualization, in a manner that allows an intuitive visual perception (e.g. a 
manner that gives a visually solid appearance), on a substantially two-dimensional field of view 
(e.g. a flat or curved plane), of a vector field (e.g. a set of topographic data representing an 
earth's surface) including a distribution of three-dimensional vectors (e.g. stereoscopic 
topographic data) that have substantially three components indicated together, or 

15 three-dimensional vectors of specified three components of multi-dimensional vectors (e.g. such 
data that have stereoscopic topographic data and geologic data together). 

The present invention further relates to a visualization processing system, a 
visualization processing method, and a visualization processing program which are adapted for 
an expression using a color tone, in place of contour lines, of altitude and inclination of 

20 concavo-convex part of terrain based on a great amount of digital image data represented by 
three-dimensional coordinates, to thereby generate a gradient-reddening stereoscopic image that 
can provide a visually solid appearance. 

Background Art 

25 For a visualization of three-dimensional vector field on a two-dimensional plane, many 

attempts have been made since ancient times. 

Most typically, there has been known a method of converting two components into 
coordinate values, and plotting an intersection thereof on a two-dimensional plane, providing 
each intersection with a note of attribute of the remaining third component (e.g. a town guide 

30 map), which however is unable to afford an easy grasp of a difference of the third component. 
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In this respect, there has been made also a graphic expression of an attributive feature 
of the third component (e.g. a town street guide map), which is still bound to the localization of 
information, and has covered successive changes in the attribute. 

To this point, there has been generally employed a method of entering a continuous 
5 feature of two components (e.g. an outline such as of a coast, river, lake, or marsh) and attribute 
isopleth lines of the third component (e.g. contour lines), which is yet difficult of the intuitive 
visible perception of an attribute variation. 

A topographic map is now supposed for more specific discussion. 

In a mesh measurement by an analyzing mapper, the terrain is lattice-like divided, 
10 having altitude values given thereto, to provide a DEM (Digital Elevation Model) data. This is 
processed in a computer for calculation of parameters concerning, e.g, the terrain's heights, 
inclination angles, slope azimuths, laplacians, valley divisions, water systems, etc., allowing for 
a distribution of calculation results over a plane to be converted into a corresponding image. 

In an airborne laser measurement, available data contain more detail information. 
15 All the data is not involved in a topographic map. 

For example, information on the height and inclilnation is extracted, to be entered as 
contour lines in the map. It however is uneasy to imagine a stereoscopic terrain therefrom. 

There is also an image provided with a stereoscopic appearance, as a hill shade lighted 
from a diagonal upside, which has an emphasized inclination in a particular direction. 
20 In this concern, there is a gray scale (tone of brightness) or a rainbow color (tone of 

hue) indicated in a terrain image, which allows an intuitive visual perception of terrain's 
geometrical features and their distribution, and is useful, but unable to give an effective visually 
solid appearance. 

Reference-1: "Japanese Patent Application Laying-Open Publication 1-46902" 
25 There is also an image processed by using either an avobeground opening or an 

underground opening as a mega filter, which allows a capture of terrain's features in a relatively 
large district, but feels something missing in visually solid appearance, particularly in local 
appearances to be visually solid. 

Reference-2: "Iwate University thesis: Indication of terrain features by an opening, the 
30 photogrammetry and remote sensing, by Ryuuzou Yokoyama, Michio Shirasawa, and Yutaka 



3 



Kikuchi (1999), vo.38, no.4, 26-34". 

Description is now made of conventional methods of providing a topographic map 
with a visually solid appearance. 

(Stereo-matching image, three-dimensional image) 
5 Basically, an image that makes use of a parallax, employing two photographs. There are 

varieties of methods, such as cases by a red/blue filter, a polarizing filter, a diffraction grating, or 
a lenticurar lenz, any of which however has to be seen in a particular direction, and needs a glass. 
Moreover, expansion as well as scaling down is difficult. 

The three-dimensional image is an image looked down in a particular direction, which is 
10 inadequate to read, as having a portion unseen if in shadow, looking small if distant, and lacking 
resolution if close. Moreover, time is necessary for image creation. 
(Indication by contour lines) 
The contour line is suitable to the indication of terrain in mountainous districts, but for 
steep inclinations (e.g. a sudden cliff part) or gentle slopes or flat lands (a plain part), the reading 
15 of topographic features takes a time due to an extreme convergence or divergence of contour 
lines having stepwise allotted heights. 

The angle of inclination as well as the orientation is to be guessed from the variation of 
spacing between contour lines. Hence, being unfit in a simple expansion or scaling, it needs a 
remake in some case. 

20 Crowded contour lines have their gaps lost, and are substituted by a legend of cliff". 

This task is complex, and constitutes an impediment to vectorization. 

Small irregularity cannot be read unless a height is given to each contour line. 

(Set of image data having two-dimensional altitude values) 
In a mapping work by aerial photographic measurement, the acquisition of 
25 information is directly made of contour lines as connected particular altitudes, having no 
altitudes given between contour lines. 

In the case of mesh measurement by an analyzing mapper or airborne laser 
measurement, the DEM data is acquired, and based thereon a two-dimensional distribution of 
contour lines is determined, whereby, although contour lines are smoothed as necessary, 
30 information else than finally contained in contour lines, e.g. information of a three-dimensional 
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geometry between contour lines, is left unused. 

This invention was made in view of the foregoing points. 

It therefore is an object of the invention to provide a visualization processing system, a 
visualization processing method, and a visualization processing program, which are adapted to 
5 visualize a vector field, with local solid attributes thereof inclusive, on a substantially 
two-dimensional field of view, in a manner that allows an intuitive visible perception. 

It also is an object of the invention to provide a visualization processing system, a 
visualization processing method, and a visualization processing program, which are adapted to 
generate a gradient reddening stereoscopic image that allows at a glance a stereoscopic grasp of 

10 terrain's heights and inclination degrees. 
Summary of Invention 

To achieve the object, a visualization processing system according to the invention is 
characterized by a first operator for mapping a vector field in a three-dimensional coordinate 
space to obtain a corresponding sequence of coordinate points, a second operator for 

15 detemrining an elevation degree in a local region of a plane connecting the sequence of 
coordinate points, a third operator for determining a depression degree in the local region of the 
plane connecting the sequence of coordinate points, a fourth operator for synthesizing the 
elevation degree and the depression degree in a weighting manner to determine an 
elevation-depression degree in the local region of the plane connecting the sequence of 

20 coordinate points, and a fifth operator for mapping the coordinate space on a two-dimensional 
plane, providing a tone indication commensurate with the elevation-depression degree to a 
region on the two-dimensional plane corresponding to a divided region of the plane connecting 
the sequence of coordinate points. 

According to the invention, a vector field is mapped in a three-dimensional coordinate 

25 space, obtaining a corresponding sequence of coordinate points, and there are determined in a 
local region of a plane connecting the sequence of coordinate points an elevation degree, that is a 
rising tendency (e.g. a ridge shaping tendency in a topographic map), and a depression degree, 
that is a sinking tendency (e.g. a valley shaping tendency in a topographic map), which are 
synthesized in an end-fit weighting manner (in a broad sense inclusive of the difference), 

30 determining in the local region an elevation-depression degree, that is a rising and sinking 
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tendency (e.g. a ridge- valley shaping tendency in a topographic map), which is tone-indicated in 
a corresponding region on a two-dimensional plane, so that the vector field can be visualized, 
with its local solid attributes inclusive, on a substantially two-dimensional plane in a manner that 
allows an intuitive visible perception. 
5 The elevation degree may preferably be defined in terms of a solid angle at one side in 

the local region of the plane connecting the sequence of coordinate points. 

The depression degree may preferably be defined in terms of a solid angle at the other 
side in the local region of the plane connecting the sequence of coordinate points. 

The visualization processing system may preferably further comprise a sixth operator 
10 for determining an inclination distribution of the plane connecting the sequence of coordinate 
points, and the fifth operator may preferably provide on the two-dimensional plane a color-toned 
indication, i.e. chroma saturation indication, of the inclination distribution (more preferably, in 
reddish colors), and for a brightness thereof, give the tone indication. 

The visualization processing system may preferably further comprise a seventh 
15 operator for connecting, among the sequence of coordinate points, those coordinate points 
equivalent of an attribute in the vector field to obtain an attribute isopleth line, and an eighth 
operator for mapping the attribute isopleth line on the two-dimensional plane given the tone 
indication. 

A visualization processing method according to the invention is characterized by a first 
20 step of mapping a vector field in a three-dimensional coordinate space to obtain a corresponding 
sequence of coordinate points, a second step of determining an elevation degree in a local region 
of a plane connecting the sequence of coordinate points, a third step of detenrdning a depression 
degree in the local region of the plane connecting the sequence of coordinate points, a fourth step 
of synthesizing the elevation degree and the depression degree in a weighting manner to 
25 determine an elevation-depression degree in the local region of the plane connecting the 
sequence of coordinate points, and a fifth step of mapping the coordinate space on a 
two-dimensional plane, providing a tone indication of the elevation-depression degree to a 
region on the two-dimensional plane corresponding to a divided region of the plane connecting 
the sequence of coordinate points. 
30 A visualization processing program according to the invention is characterized in that 
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the program is functionable to have a computer execute a first process for mapping a vector field 
in a three-dimensional coordinate space to obtain a corresponding sequence of coordinate points, 
a second process for determining an elevation degree in a local region of a plane connecting the 
sequence of coordinate points, a third process for determining a depression degree in the local 
5 region of the plane connecting the sequence of coordinate points, a fourth process for 
synthesizing the elevation degree and the depression degree in a weighting manner to determine 
an elevation-depression degree in the local region of the plane connecting the sequence of 
coordinate points, and a fifth process for mapping the coordinate space on a two-dimensional 
plane, providing a tone indication of the elevation-depression degree to a region on the 
10 two-dimensional plane corresponding to a divided region of the plane connecting the sequence 
of coordinate points. 

According to the invention, a visualization processing system for generating a gradient 
reddening stereoscopic image is characterized by a database having stored therein a multiplicity 
of digital data provided with three-dimensional coordinates, and a computer which comprises a 

15 means for generating a stereoscopic contour image having contour lines connecting 
three-dimensional coordinates of digital data having identical Z values, a means for meshing 
intervals between contour lines, a means for allocating focused points to meshes, determining an 
average of differences in Z value between a respective mesh given a focused point and 
neighboring meshes, a means for generating a gradient reddening image having assigned to the 

20 mesh given the focused point a red tone commensurate with a degree in magnitude of a 
difference of the average, a means for generating a gray scale image having a varied brightness 
depending on a ridge- valley shaping tendency of the mesh gi ven the focused point, and a means 
for performing a multiplying synthesis of the gradient reddening image and the gray scale image, 
to display on a screen a gradient reddening stereoscopic image representing degrees of gradient 

25 and degrees of height in color. 

According to the invention, a visualization processing method for generating a gradient 
reddening stereoscopic image is characterized by a step of generating a stereoscopic contour 
image having contour lines connecting three-dimensional coordinates of digital data having 
identical Z values, a step of meshing intervals between contour lines, a step of allocating focused 

30 points to meshes, determining an average of differences in Z value between a respective mesh 
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given a focused point and neighboring meshes, a step of generating a gradient reddening image 
having assigned to the mesh given the focused point a red tone commensurate with a degree in 
magnitude of a difference of the average, a step of generating a gray scale image having a varied 
brightness depending on a ridge-valley shaping tendency of the mesh given the focused point, 
5 and a step of performing a multiplying synthesis of the gradient reddening image and the gray 
scale image, to display on a screen a gradient reddening stereoscopic image representing degrees 
of gradient and degrees of height in color. 

According to the invention, a visualization processing program for generating a 
gradient reddening stereoscopic image is characterized in that the program is adapted to have a 

10 computer function as a means for reading a multiplicity of digital data provided with 
three-dimensional coordinates, a means for generating a stereoscopic contour image having 
contour lines connecting three-dimensional coordinates of digital data having identical Z values, 
a means for meshing intervals between contour lines, a means for allocating focused points to 
meshes, determining an average of differences in Z value between a respective mesh given a 

15 focused point and neighboring meshes, a means for generating a gradient reddening image 
having assigned to the mesh given the focused point a red tone commensurate with a degree in 
magnitude of a difference of the average, a means for generating a gray scale image having a 
varied brightness depending on a ridge-valley shaping tendency of the mesh given the focused 
point, and a means for performing a multiplying synthesis of the gradient reddening image and 

20 the gray scale image, to display on a screen a gradient reddening stereoscopic image 
representing degrees of gradient and degrees of height in color. 

Brief Description of the Drawings 

The above and further features, functions, and effects of the invention become more 
25 apparent from the following description of best modes for carrying out the invention, when the 
same is read in conjunction with the accompanying drawings, in which: 

Fig. 1 is a block diagram of a visualization processing system according to a first mode 
of embodiment of the invention; 

Fig. 2 is a flowchart showing a processing procedure and process results of the 
30 visualization processing system of Fig. 1 ; 
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Fig. 3 is a detail of part A of the flowchart of Fig. 2; 
Fig. 4 is a detail of part B of the flowchart of Fig. 2; 
Fig. 5 is a detail of part C of the flowchart of Fig. 2; 
Fig. 6 is a detail of part D of the flowchart of Fig. 2; 
5 Fig. 7 is a detail of part E of the flowchart of Fig. 2; 

Fig. 8 is a detail of part F of the flowchart of Fig. 2; 
Fig. 9 is a detail of part G of the flowchart of Fig. 2; 
Fig. 10 is a detail of part H of the flowchart of Fig. 2; 

Fig. 11 is a block diagram of a visualization processing system for generating a 
10 gradient reddening stereoscopic image according to a second mode of embodiment of the 
invention; 

Fig. 12 is an illustration of a laser measurement in the visualization processing system 
of Fig. 11; 

Fig. 13 is an illustration of an eight-directional array in the visualization processing 
15 system of Fig. 1 1 ; 

Fig. 14 is an illustration of principles of an aboveground opening and an underground 
opening in the visualization processing system of Fig. 1 1 ; 

Fig. 15 is a diagram illustrating principal patterns of aboveground opening and 
underground opening in the visualization processing system of Fig. 1 1 ; 
20 Fig. 16 is a stereoscopic illustration of an aboveground opening and an underground 

opening in the visualization processing system of Fig. 1 1 ; 

Fig. 17 is a diagram illustrating a sample point and a distance of aboveground opening 
and underground opening in the visualization processing system of Fig. 1 1 ; 

Fig. 18 is a diagram illustrating gray scale allocation in the visualization processing 
25 system of Fig. 1 1 ; 

Fig. 19 is a block diagram of an convexity -emphasized image creator, a 
concavity-emphasized image creator, and a first synthesizer of the visualization processing 
system of Fig. 11; 

Fig. 20 is a block diagram of an inclination-emphasized image creator, and a second 
30 synthesizer of the visualization processing system of Fig. 1 1 ; 
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Fig. 21 is an illustration of a generation step of a synthetic image of an aboveground 
opening image and an underground opening image in the visualization processing system of Fig. 

ll; 

Fig. 22 is an illustration of steps for generation of a gradient reddening stereoscopic 
5 image in the visualization processing system of Fig. 1 1 ; 

Fig. 23 is a view of a stereoscopically visualized image of Aokigahara, Mt. Fuji, 
obtained by the visualization processing system of Fig. 1 1 ; 

Fig. 24 is a contour-lined topographic map based on an aerial mapping measurement 
of a southern region of a Tenjiny ama ski area, Aokigahara, Mt. Fuji; 
10 Fig. 25 is a contour-lined topographic map based on laser measurement data of the 

same region; 

Fig. 26 is a magnified view of a stereoscopically visualized image of the same region; 

Fig. 27 is a magnified view of a stereoscopically visualized image of another region; 

Fig. 28 is a magnified view of a stereoscopically visualized image of another region; 
15 Fig. 29 is an X-ray fluoroscopic view of a human body; 

Fig. 30 is a graduation-tinted slice image of the fluoroscopic view of Fig. 29; 

Fig. 31 is a red-toned elevation-depression degree distribution image obtained by 
processing the fluoroscopic view of Fig. 29 by the visualization processing system of Fig. 1 ; and 

Fig. 32 is a synthetic image of the image of Fig. 30 and the image of Fig. 31 
20 superposed thereon. 

Best Modes for Carrying Out the Invention 

There will be described below best modes for carrying out the invention. 

25 (First mode of embodiment) 

First, description is made of a first mode of embodiment of the invention, with 
reference to Fig. 1 to Fig. 10. Fig. 1 is a block diagram of a visualization processing system 
VPS1 according to this mode of embodiment, Fig. 2, a flowchart showing process procedures 
PI to P8 and process results A to I of the system VPS1, and Fig. 3 to Fig. 10, details of principal 

30 process results A to H, respectively. 
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As shown in Fig. 1, the visualization processing system VPS1 has a central 
information processing unit (CPU) 51 composed as an adequate combination of a workstation, a 
processor, or a microcomputer, and logics, registers, etc., an information input section 54 
including a keyboard (KB), a mouse, interactive software switches, external communication 
5 channels, etc., for inputting necessary control and/or operational information to the central 
information processing unit 51, an information output section 55 including a display, a printer, 
external communication channels, etc., for performing indication and/or transmission, in a wide 
sense, of information output from the central information processing unit 51, a first memory 52 
such as a read-only memory (ROM) having stored therein an operating system, application 

10 programs, and the like to be read to the central information processing unit 51, and a second 
memory 53 such as a random access memory (RAM) for storing information to be occasionally 
processed in the central information processing unit 51 as well as information occasionally 
written from the central information processing unit 51. The first and second memories 52, 53 
may preferably be integrated or subdivided in a suitable manner. 

15 The first memory 52 has stored therein a visualization processing program 60 operable 

on a prescribed application. In this mode of embodiment, the visualization processing program 
60 is composed of a first to an eighth processing file 61 to 68 including program groups to be 
read in the central information processing unit 51, where they are adapted to execute eight 
visualization processes PI to P8 shown in Fig. 2, respectively, while the classification and 

20 allocation of those program groups can be set in a voluntary manner. 

In the second memory 53 is stored a vector field 70 constituting an object of process of 
the visual processing program 60. The vector field 70 may be a finite set of (a total number of N) 
information vectors having one or more components allowing extraction of substantially three or 
more kinds of information. In this mode of embodiment, each vector is given as a 

25 two-component vector that contains, for a focused point representing a minute finite-division 
region of a ground surface of Mt. Fuji, an identification (Id) number allowing confirmation to be 
made of information on the latitude and information on the longitude by a reference table, and a 
height difference relative to a neighboring focused point or reference point of triangulation. 

The first processing file 61 is adapted to calculate, from an identification number Id n 

30 and a height difference of a two-component vector V n processed as an n-th (n = 1 to N), the 
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longitude x n , latitude y n , and sea level altitude Zn, and associate their values with a corresponding 
coordinate point Qn = {X n = x n , Y n = y n , Zn = Zn} in a virtual three-dimensional (3D) X-Y-Z 
orthogonal coordinate space 80 stored in the second memory 53, i.e., store the identification 
number Id of the vector V n in a storage region in the memory 53 correspondent to the coordinate 
5 point Qn, to thereby map the vector V n into the coordinate space 80. This is made for the total 
number of N vectors, whereby the vector field 70 is mapped in the coordinate space 80 (Fig. 2, 
process PI). 

The first processing file 61 is further adapted to determine by the method of least 
squares or the like a curved plane S connecting with a necessary smoothness a sequence of a 
10 total number of N or an adequate smaller number of Id-numbered coordinate points { Qn* n < N } 
within the coordinate space 80, divide it into a total number of M {M<N} minute plane 
regions {S m : m<M], defining their focused points Qm, and store relevant information in the 
memory 53. 

A second processing file 62 is adapted to verify, for a respective plane region S m , a 
15 local region Lj n + at an obverse side (Z+ side) of the curved plane S residing within a prescribed 
radius from a focused point Qm thereof, and determine a degree of openness defined thereby (i.e. 
a see-through solid angle to the heaven end or a second-order differential value equivalent 
thereto) m + about the focused point Qm (Fig. 2, process P2), storing it as an elevation degree 
of the plane region S m . Fig. 3 shows, as a process result A, an image in which the elevation 
20 degree *F m + is tone-indicated over an entirety of the curved plane S. This image A clearly 
indicates a ridge side of terrain, i.e., a convexity (of the curved plane S) like an evident 
convexity. 

A third processing file 63 is adapted to verify, for the plane region S m , a local region 
Lm" at a reverse side (Z- side) of the curved plane S residing within the prescribed radius from 

25 the focused point Qm, and determine a degree of openness defined thereby (i.e. a see-through 
solid angle to the earth end or a second-order differential value equivalent thereto) ¥ m about 
the focused point Qm (Fig. 2, process P3), storing it as a depression degree of the plane region S m . 
Fig. 5 shows, as a process result C, an image in which the depression degree *F m ~ is 
tone-indicated over an entirety of the curved plane S. This image C clearly indicates a valley side 

30 of terrain, i.e., a concavity (of the curved plane S) like an evident concavity. It should be noted 
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that this image C does not constitute a simple reverse of the image A. 

A fourth processing file 64 is adapted to synthesize, for the plane region S m , the 
elevation degree ¥ m + and the depression degree ¥ m in a weighting manner (w +v F m + + 
w" *F m ") with a sharing proportion w + : w" (w + + w" = 0) determined in an end-fit manner (that is, 
5 depending on which of ridge and valley is to be put above), thereby deteninining a stereoscopic 
effect to be brought about the focused point Q m by a local region Lm (L m + , Lm") at obverse and 
reverse of the curved plane S residing within the prescribed radius (Fig. 2, process P4), storing it 
as an elevation-depression degree ¥ m of the plane region S m . Fig. 4 shows, as a process result 
B, an image in which the elevation-depression degree *F m is tone-indicated over an entirety of 

10 the curved plane S. This image B clearly indicates a convexity (of the curved plane S) like an 
evident convexity, and a concavity like an evident concavity, thereby defining ridge and valley 
of terrain, with an intensified visually solid feeling. It is noted, for the image B, the weighting in 
synthesis is w + = -w = 1. 

Description is now made of a sixth processing file 66. This file 66 is adapted to 

15 determine, for the plane region S m , a maximum degree of inclination G m thereof (or a first-order 
differential value equivalent thereto) directly or indirectly via the method of least squares (Fig. 2, 
process P6), storing it as an inclination G m of the plane region S m . Fig. 6 shows, as a process 
result D, (an achromatic indication image of) an image in which the inclination G m is 
color-toned to be indicated in a red-spectral color R over an entirety of the curved plane S. This 

20 image D also has the effect of visually projecting a stereoscopic appearance of terrain (that is the 
curved plane S). 

A fifth processing file 65 is adapted for mapping (Fig. 2, process P5) the 
three-dimensional coordinate space 80, together with relevant information thereof (*F m , G m , R), 
onto a two-dimensional plane 90 in the information output section 55, to thereby provide, to a 

25 region 90 m on the two-dimensional plane 90 corresponding to the region S m as a division of the 
plane S connecting a sequence of coordinate points Qm, an R color-toned indication of the 
inclination G m , and for a brightness of the R color tone, a tone indication commensurate with the 
elevation-depression degree *F m . Fig. 8 shows, as a process result F, (an achromatic indication 
image of) this image. In this image F, terrain (that is the curved plane S) has a visually solid 

30 appearance. 
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An image E of Fig. 7 shows a result of a mapping (process P5), by the processing file 
65 onto the two-dimensional plane 90, of information on the image D (i.e.R color tone 
indication of inclination GJ and information on elevation-depression degrees (i.e. elevation 
degrees *F m + ) corresponding to the image A, where ridge part is emphasized. 
5 An image G of Fig. 9 shows a result of a mapping (process P5), by the processing file 

65 onto the two-dimensional plane 90, of information on the image D (Le. R color tone 
indication of inclination Gm) and information on elevation-depression degrees (i.e. depression 
degrees *P m ) corresponding to the image C, where valley part is emphasized. 

A seventh processing file 67 is adapted to determine attribute isopleth lines (in this 

10 mode of embodiment, terrain's isometric contour lines and shape contour lines) Ea connecting, 
among the sequence of coordinate points Qn, those coordinate points Qn equivalent of an 
attribute (in this mode of embodiment, sea level altitude Zn) extracted from components of 
vectors V n of the vector 70 field, storing them, and to output or indicate, as necessary (Fig. 2, 
process P7). Fig. 2 shows an indication process result I of the same. This result I also contributes 

15 to the grasp of a stereoscopic configuration of terrain (that is the curved plane S). 

The eighth processing file 68 is adapted to map or output for display the 
three-dimensional space 80, together with relevant information ( *P m , G m , R) thereof, onto the 
two-dimensional plane 90, mapping or outputting for display the attribute isopleth lines Ea (Fig. 
2, process P8). Fig. 10 shows, as a process result H, (an achromatic indication image of) a 

20 display image of the same. In this image H also, terrain (that is the curved plane S) has a visually 
solid appearance. 

Therefore, the visualization processing system VPS1 according to this mode of 
embodiment comprises a first operator (61) for mapping a vector field 70 in a three-dimensional 
coordinate space 80 to obtain a corresponding sequence of coordinate points Qm, a second 

25 operator (62) for determining an elevation degree *F m + in a local region Lm + of a plane S 
connecting the sequence of coordinate points, a third operator (63) for determining a depression 
degree ¥ m in a local region of the plane S connecting the sequence of coordinate points, a 
fourth operator (64) for synthesizing the elevation degree and the depression degree in a 
weighting manner to determine an elevation-depression degree ¥ m in a local region of the 

30 plane S connecting the sequence of coordinate points, and a fifth operator (65) for mapping the 



14 

coordinate space 80 on a two-dimensional plane 90, providing a tone indication commensurate 
with the elevation-depression degree to a region 90 m on the two-dimensional plane 90 
corresponding to a divided region S m of the plane S connecting the sequence of coordinate 
points. The operator means herein an element, a set thereof, or a means for executing an 
5 operation process command or operation process function programmed or set in advance. 

The visualization processing system VPS1 further comprises a sixth operator (66) for 
determining an inclination G m distribution of the plane S connecting the sequence of coordinate 
points, and the fifth operator (65) is adapted to provide on the two-dimensional plane 90 a 
color-toned indication of the inclination distribution, in a red spectral color R, and for a 

10 brightness thereof, give the tone indication. 

The visualization processing system VPS1 further comprises a seventh operator (67) 
for connecting, among the sequence of coordinate points, those coordinate points equivalent of 
an attribute in the vector 70 field to obtain attribute isopleth lines Ea, and an eighth operator (68) 
for mapping the attribute isopleth lines Ea on the two-dimensional plane 90 given the tone 

15 indication. 

According to this mode of embodiment, a vector field 70 can be visualized on a 
substantially two-dimensional plane 90 in a manner that allows an intuitive visual perception, 
with a local stereoscopic attribute thereof inclusive. 

A visualization processing method according to this mode of embodiment comprises a 

20 first step PI of mapping a vector field in a three-dimensional coordinate space to obtain a 
corresponding sequence of coordinate points, a second step P2 of determining an elevation 
degree in a local region of a plane connecting the sequence of coordinate points, a third step P3 
of determining a depression degree in the local region of the plane connecting the sequence of 
coordinate points, a fourth step P4 of synthesizing the elevation degree and the depression 

25 degree in a weighting manner to determine an elevation-depression degree in the local region of 
the plane connecting the sequence of coordinate points, and a fifth step P5 of mapping the 
coordinate space on a two-dimensional plane, providing a tone indication of the 
elevation-depression degree to a region on the two-dimensional plane corresponding to a 
divided region of the plane connecting the sequence of coordinate points. 

30 A visualization processing program 60 according to this mode of embodiment is 
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functionable to have a central information processing unit 51 execute a first process PI for 
mapping a vector field in a three-dimensional coordinate space to obtain a corresponding 
sequence of coordinate points, a second process P2 for determining an elevation degree in a 
local region of a plane connecting the sequence of coordinate points, a third process P3 for 
5 determining a depression degree in the local region of the plane connecting the sequence of 
coordinate points, a fourth process P4 for synthesizing the elevation degree and the depression 
degree in a weighting manner to determine an elevation-depression degree in the local region of 
the plane connecting the sequence of coordinate points, and a fifth process P5 for mapping the 
coordinate space on a two-dimensional plane, providing a tone indication of the 
10 elevation-depression degree to a region on the two-dimensional plane corresponding to a 
divided region of the plane connecting the sequence of coordinate points. 

(Second mode of embodiment) 

Next, description is made of a second mode of embodiment of the invention, with 

15 reference to Fig. 1 1 to Fig. 28. 

Fig. 11 is a block diagram of a visualization processing system VPS2 including a 
gradient reddening stereoscopic image generator 4 according to this mode of embodiment, Fig. 
12, an illustration of a laser measurement in the system VPS2, Fig. 13, an illustration of an 
eight-directional array, Fig. 14, an illustration of principles of an aboveground opening and an 

20 underground opening, Fig. 15, a diagram illustrating principal patterns of aboveground opening 
and underground opening, Fig. 16, a stereoscopic illustration of an aboveground opening and an 
underground opening, Fig. 17, a diagram illustrating a sample point and a distance of 
aboveground opening and underground opening, Fig. 18, a diagram illustrating gray scale 
allocation, Fig. 19, a block diagram of an convexity-emphasized image creator, a 

25 concavity-emphasized image creator, and a first synthesizer, Fig. 20, a block diagram of an 
inclination-emphasized image creator, and a second synthesizer, Fig. 21, an illustration of a 
generation step of a synthetic image of an aboveground opening image and an underground 
opening image, and Fig. 22, an illustration of steps for generation of a gradient reddening 
stereoscopic image. 

30 Further, Fig. 23 is a view of a stereoscopically visualized image of Aokigahara, Mt. 
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Fuji, obtained by the visualization processing system VPS2, Fig. 24, a contour-lined topographic 
map based on an aerial mapping measurement of a southern region of a Tenjinyama ski area, 
Aokigahara, Mt. Fuji, Fig. 25, a contour-lined topographic map based on laser measurement 
data of the same region, Fig. 26, a magnified view of a stereoscopically visualized image of the 
5 same region, Fig. 27, a magnified view of a stereoscopically visualized image of another region, 
and Fig. 28, a magnified view of a stereoscopically visualized image of another region. 

This mode of embodiment is configured to determine an inclination corresponding to 
the inclination G m in the first mode of embodiment, an aboveground opening corresponding to 
the elevation degree *F m + in the first mode of embodiment, and an underground opening 
10 corresponding to the depression degree ^ m in the first mode of embodiment, as three 
parameters based on DEM (Digital Elevation Model) data, and store their flat plane distribution 
as a gray scale image. 

A difference image between aboveground and underground is entered in a gray, and 
the inclination, in a red channel, for creation of a pseudo-color image, to thereby effect indication 

15 to be whitish at a ridge or crest part, and blackish in valley or howe, while it gets red as inclined 
part becomes steep. Such combination of indication allows creation of an image (hereafter 
sometimes called a stereoscopic reddened map) to be solid in appearance even on a single sheet. 

In other words, for a map to be solid, this mode of embodiment employs a stereoscopic 
indication method, in which intervals between contour lines are meshed, and for a respective 

20 mesh, a difference to a neighboring mesh, that is, an inclination is indicated by a red color tone, 
and whether high or low in comparison with periphery is indicated by a gray scale. This 
corresponds to the elevation-depression degree *F m in the first mode of embodiment, which is 
referred to as a ridge-valley shaping tendency in this mode of embodiment, suggesting a brighter 
place to be (ridge-like) higher than periphery, and a darker place to be (valley-like) lower than 

25 periphery, of which brightness and darkness are multiplication-synthesized to thereby generate a 
solid appearance. 

Fig. 11 shows a schematic arrangement of the gradient reddening stereoscopic image 
generator 4 according to this mode of embodiment. As shown in Fig. 11, the gradient reddening 
stereoscopic image generator 4 of this mode of embodiment is provided with a computer 
30 function described below. The gradient reddening stereoscopic image generator 4 further has a 
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variety of databases connected thereto. 

A database 1 has a laser data Ri stored therein. For the laser data (Rx, Ry, Rz: 
coordinates of laser data referred to by addition of R), as illustrated in Fig. 12, by an aircraft 
flying level in the sky above a target district (preferably, a digital camera pickup range), a laser 
5 beam is emitted downward, and from the time required to go and back, aircraft's location and 
posture, and angle of emission, ground surface's x, y, z are determined by (computer) calculation 
and stored. Employed for grasp of the aircraft's location is a GPS (not shown), and for grasp of 
the posture, an IMU. 

The laser emitter (unshown at Z) is adapted for 33,000 shots per second, allowing 
10 acquisition of altitude points (Rx, Ry, Rz) in a density of one point every 80cm. 

In the case a plurality of reflected pulses are measured for one shot of laser emission, 
data of a final reflection is employed and stored. 

Further, by checking received laser data for a tendency of distribution, points spiked 
higher than periphery are removed as being laser data of trees having failed to be passed, in 
15 addition to removal of laser data else than tree, such as a house, automobiles, or bridge. 
Therefore, the database 1 simply stores a laser data Ri of ground surface. 

A database 2 has stored therein at least a contour map Hi (1/25,000: with numbered 
contour lines) of the digital camera pickup range. Further, the contour map is provided with 
characteristic point coordinates (Hx, Hy, Hz; contour map data). 
20 Further, a database 3 has stored therein a stereo-matching data Mi. The 

stereo-matching data Mi is created as a stereoscopic image from two air photos that have picked 
up an identical area. For example, a known building is extracted from the two photos, and a side 
of the building is given a Z value for solidization (Mx, My, Mz), which constitutes a reference 
for Z values to be given to others. 
25 A DEM data creator 6 reads laser data Ri of the database 1 , generating a contour map 

in which respective altitudes identical in value are connected, and creates a TIN for the contour 
map, to restore a ground surface. It then determines heights of crossing points between the TIN 
and respective lattice points, to produce DEM (Digital Elevation Model) data. 

Further, the DEM data creator 6 reads a contour map Hi stored in the database 2, and 
30 generates the TIN mutually connecting contour lines, which is converted to the above-noted 
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DEM data. 

Next, description is made of a DEM data to be employed in this mode of embodiment. 
For example, for a "numeric map 50-m mesh (altitude)", from meshes of a 1/25,000 topographic 
map vertically and horizontally divided into 200 equal parts (by mesh pitches of 2.25 seconds in 
5 latitude-line direction and 1.50 seconds in longitude-line direction), their central altitudes are 
read at intervals of 1 m, to be arrayed two-dimensional. 

Further, in this mode of embodiment, the gradient reddening stereoscopic image 
creator 4 includes an aboveground opening data creator 9, an underground opening data creator 
10, a gradient calculator 8, a convexity-emphasized image creator 11, a concavity-emphasized 
10 image creator 12, an inclination emphasizer 13, a first synthesizer 14, and a second synthesizer 
15. 

In this mode of embodiment, the concept of an opening is used. The opening is a 
quantified degree by which a spot in concern is convex aboveground or concave underground in 
comparison with surroundings. In other words, as illustrated in Fig. 14, an aboveground opening 

15 is defined as an extent of the sky to be seen within a range of a distance L from a focused sample 
point, and an underground opening is defined as an extent under the ground, within a range of 
the distance L, when taking a survey in the soil in a handstand position. 

The openings depend on the distance L and a surrounding terrain. Fig. 15 illustrates, 
for 9 kinds of principal terrains, their aboveground opening and underground opening by graphs 

20 of octagons representing an aboveground angle and an underground angle in respective 
azimuths. Generally, the aboveground opening increases as the spot is projected higher from the 
surrounding, and has a value to be large at a crest or ridge, but small in a howe or at a valley 
bottom. To the contrary, the underground opening increases as the spot is cut underground lower, 
and has a value to be large in a howe or at a valley bottom, but small at a crest or ridge. Actually, 

25 even in the range of distance L, a variety of principal terrains are weaved, so that the octagonal 
graphs of aboveground angle and underground angle are frequendy deformed, giving varieties 
of opening values. 

As will be described, D <j> L and D y/ L have non-increasing characteristics to L, and 
<p L and *F L have non-increasing characteristics to L, accordingly. 
30 Further, the opening diagram permits extraction of information to fit to the terrain scale, 
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by specification of a calculated distance, allowing for an indication free from dependency on 
directional and local noises. 

That is, excellent extraction of ridge line and valley line allows abundant geographical 
and geological information to be read: as illustrated in Fig. 16, on a DEM data (ground surface: 
5 solid: Fig. 16(a)) within a fixed range, a spot A in concern is set, which is connected with a point 
B to be a greatest peak when viewed in any one of eight directions therefrom, by a straight line 
LI that defines, to a horizontal line, an angular vector 0 i to be determined. In this way, over the 
eight directions, angular vectors are determined, of which an average is called an aboveground 
opening # i; and on an invert DEM data (Fig. 16(c)), which is a reverse of a solid (Fig. 16(b)) 

10 having an air layer pushed on the DEM data of the fixed range (ground surface: solid), the spot 
A in concern is connected with a point C (corresponding to a deepest place) to be a greatest peak 
when viewed in any one of eight directions therefrom, by a straight line L2 that defines, to the 
horizontal line, an angle to be determined. This angle is determined over the eight angles, and 
averaged, which is called an underground opening y/ i. 

15 That is, the aboveground opening data creator 9 produces, on a DEM data included 

within a range of a fixed distance from a focused point, a terrain section in each of eight 
directions, and determines a maximum value (in view of a plumb direction) among inclinations 
of connection lines (LI of Fig. 16(a)) between the focused point and respective terrain points. A 
process like this is executed in the eight directions. The angle of inclination is an angle from the 

20 zenith (90 degrees on a flat, above 90 degrees at a ridge or crest, or below 90 degrees at a valley 
bottom or in a howe); and the underground opening data creator 10 produces, within a range of 
the fixed distance from the focused point on an invert DEM data, a terrain section in each of 
eight directions, and determines a maximum value among inclinations of connection lines 
between the focused point and respective terrain points (to be a minimum value when looking 

25 L2 in the plumb direction in a three-dimensional figure of a ground surface of Fig. 16(a)). A 
process like this is executed in the eight directions. 

When looking L2 in the plumb direction in the three-dimensional figure of the ground 
surface of Fig. 16(a), the angle y/ i has 90 degrees on a flat, under 90 degrees at a ridge or crest, 
or above 90 degrees at a valley bottom or in a howe. 

30 In other words, for the aboveground opening and underground opening, as illustrated 
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in Fig. 17, two base points A (iA, jA, HA) and B (iB, jB, HB) are supposed. As the sample 
spacing is 1 m, the distance between A and B is given, such that 

Fig. 17(a) illustrates a relationship between sample points A and B, relative to an 
5 altitude 0 m as a reference. An elevation angle 0 the sample point A has to the sample point B 
is given, such that 

0 = tan~ ] {(HB- HA)/ P. 
The sign of 0 is positive for ® HA<HB, or negative for (2) HA>HB. 

A set of sample points residing in an azimuth D within a range of a distance L from a 
10 focused sample point is denoted DSL, which will be called "a D-L set of a focused sample 
point". Letting now 

D J3 L: a maximum value among elevation angles for respective elements of DSL of a 
focused sample point, and 

D 5 L: a minimum value among elevation angles for respective elements of DSL of a 
15 focused sample point (refer to Fig. 17(b)), 
the following definitions are given. 

Definition 1: an aboveground angle and an underground angle of a D-L set of a 
focused sample point shall mean respectively such that 

D ( 2)L = 90-D > 3L,and 
20 D^L = 90 + D£L. 

D (j> L means a maximum value of a zenith angle in which the sky in an azimuth D can 
be seen within a distance L from a focused sample point. A generally called horizon angle 
corresponds to the aboveground angle as L is an infinity. And, D y/ L means a maximum value 
of a nadir angle in which the soil in an azimuth D can be seen within a distance L from a focused 
25 sample point. As L is increased, the number of sample points belonging to DSL increases, 
whereto D f3 L has a non-decreasing characteristic, and on the contrary, D 8 L has a 
non-increasing characteristic. Therefore, D (f> L as well as D^l has a non-increasing 
characteristic to L. 

In the geodesy, the high angle is a concept defined to a horizontal plane as a reference 
30 passing a focused sample point, and strictly, not coincident with 6. Further, for a strict 
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discussion of aboveground angle and underground angle, the curvature of the earth should also 
be considered, and the definition 1 is not always an exact description. The definition 1 is a 
concept defined for a geomorphic analysis assumed to be made by using DEM to the last. 

The aboveground angle and the underground angle have been concepts for a specified 
5 azimuth D, which will be expanded by introducing the following definition. 

Definition II: An aboveground opening and an underground opening of a distance L 
from a focused sample point shall mean respectively such that 
OL=(0^L + 45^L + 90^L+135^L+180^L 

+ 225 <f> L + 270 <j) L + 3 1 5 <f) L) / 8, and 
10 ¥L=(0^L + 45^L + 90^L+ 135 i//L+ 180^ L 

+ 225^L + 270^L+315^L)/8. 
The aboveground opening represents an extent of the sky to be seen within a range of 
the distance L from the focused sample point, and the underground opening represents an extent 
under the ground, within a range of the distance L, when taking a survey in the soil in a 
15 handstand position (refer to Fig. 14). 

The gradient calculator 8 is adapted to mesh a DEM data into squares, and determine 
an average gradient of surfaces of squares neighboring a focused point on the meshes. The 
neighboring squares are four in number, any one of which is chosen as a focused square. Then, 
an average inclination and heights at four comers of the focused square are determined. The 
20 average inclination is a gradient of surface approximated from four points by using the method 
of least squares. 

The convexity-emphasized image creator 1 1 has a first gray scale for indicating a ridge 
and a valley bottom by brightness, as illustrated in Fig. 18(a), and is adapted, every time when 
the aboveground opening data creator 9 determines an aboveground opening (as an average 
25 angle when a range of L is seen in eight directions from a focused point: an index for judgment 
of whether the residing place is high), to calculate a brightness (luminance) commensurate with 
the aboveground opening 0 \. 

For example, for values of aboveground opening falling within a range of about 40 
degrees to 120 degrees, the first gray scale is associated with a range of 50 degrees to 110 
30 degrees, which is allotted to 255 tones. 
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That is, toward a ridge part (convex), the place has a greater value of aboveground 
opening, and becomes white in color. 

Then, as shown in Fig. 19, at the convexity-emphasized image creator 11, a convexity 
emphasizing color allotment process 20 reads an aboveground opening image data Da, and 
5 allots, to a mesh region having a focused point (coordinate) (in the case contour lines connecting 
identical Z values of DEM data are meshed (e.g. 1 m) into squares, and a focused point is set to a 
point at any of four comers of mesh) a color data based on the first gray scale, which is stored (as 
an aboveground opening image data Dpa) in an aboveground opening file 21. Then, a tone 
corrector 22 stores, in a file 23, the aboveground opening image data Dpa having color tones 
10 thereof inverted, as an aboveground opening layer Dp. That is, there is obtained an aboveground 
opening layer Dp adjusted for the ridge to be indicated white. 

A concavity extractor 12 has a second gray scale for indicating a valley bottom and a 
ridge by brightness, as illustrated in Fig. 18(b), and is adapted, every time when the underground 
opening data creator 10 determines an underground opening y/i (as an average in eight 
15 directions from a focused point), to calculate a brightness commensurate with the underground 
opening y/i. 

For example, for values of underground opening falling within a range of about 40 
degrees to 120 degrees, the second gray scale is associated with a range of 50 degrees to 110 
degrees, which is allotted to 255 tones. 
20 That is, toward a valley bottom part (concave), the place has a greater value of 

underground opening, and becomes black in color. 

Then, as shown in Fig. 19, at the concavity-emphasized image creator 12, a concavity 
emphasizing color allotment process 25 reads an underground opening image data Db, and 
allots, to a mesh region having a focused point (coordinate) (in the case contour lines connecting 
25 identical Z values of DEM data are meshed (e.g. 1 m) into squares, and a focused point is set to a 
point at any of four comers of mesh) a color data based on the second gray scale, which is stored 
in an underground opening file 26. Then, a tone correction process 27 makes a correction of 
color tones of the underground opening image data Db. 

If the color is toned excessively black, it is set to a tone of color according to a 
30 corrected tone curve. This is called an underground opening layer Dq, and stored in a file 28. 
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The inclination-emphasized image creator 13 has a third gray scale for indicating a 
degree of inclination by brightness, as illustrated in Fig. 18(c), and is adapted, every time when 
the gradient calculator 8 determines an inclination degree (as an average in four directions from 
a focused point), to calculate a brightness (luminance) commensurate with the inclination 
5 degree. 

For example, for values of inclination a i falling within a range of about 0 degree to 
70 degrees, the third gray scale is associated with a range of 0 degree to 50 degrees, which is 
allotted to 255 tones. That is, 0 degree to be white, and 50 degrees to be black. The color is 
blackened as the spot has a greater inclination a . 

10 Then, as shown in Fig. 20, at the inclination-emphasized image creator 13, an 

inclination emphasizing color allotment process 30 stores, in a file 31, a difference image 
between the underground opening image data Db and the aboveground opening image data Da, 
as a gradient image Dra. 

At this time, a color data based on the third gray scale is allotted to a mesh region 

15 having a focused point (coordinate) (in the case contour lines connecting identical Z values of 
DEM data are meshed (e.g. 1 m) into squares, and a focused point is set to a point at any of four 
corners of mesh). Then, a reddening process has an RGB color mode function for emphasizing 
R. That is, there is obtained in a file 33 a gradient-emphasized image Dr that has an emphasized 
red as the inclination is greater. 

20 The first synthesizer 14 is adapted for a multiplication of the aboveground opening 

layer Dp and the underground opening layer Dq to obtain a synthetic image Dh (Dh = Dp + Dl) 
thereby synthesized. At this time, a balance of both is adjusted to avoid collapsing valley part. 



The above-noted "multiplication" is a term in a layer mode on a photoshop, that 
corresponds to an OR operation for numeric process. 



Brightness 


Gray scale 


i-1 


i-2 


i-3 


Lighter as higher 


aboveground 


36 


52 


45 


Darker as higher 


underground 


32 


48 


61 




Total 


68 


100 


106 



25 In the balance adjustment, for allocation of value between aboveground and 

underground at a spot, the ground surface is cut out by a fixed radius (L/2) about the spot. 
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Assuming an entirety of sky to be uniform in brightness, the extent of sky looked up 
from a surface of ground gives a brightness of the ground surface. 

That is, an aboveground opening constitutes the brightness. However, assuming light 
spiking around, a value of underground opening should be considered. 
5 Depending on how the ratio of both is settled, an arbitrary modification can be 

achieved, such as for an emphasis at a ridge part of terrain, for example. For a desirable 
emphasis of a terrain in a valley, the value of b may be increased. 
Index of brightness = a X aboveground opening 
-b X underground opening, 

10 where a + b=l. 

That is, as illustrated in Fig. 21, the aboveground opening layer Dp (with ridge 
emphasized white) and the underground opening layer Dq (with bottom emphasized black) are 
multiplication-synthesized to obtain a synthetic image gray-toned for indication (Dh = Dp + 
Dl). 

15 On the other hand, the second synthesizer is adapted for synthesizing the 

gradient-emphasized image Dr in the file and the synthetic image Dh obtained by a synthesis at 
the first synthesizer, to have a stereoscopic reddening image Ki emphasized at the ridge in red 
color, for indication on the display. 

That is, as illustrated in Fig. 22, there is obtained the gradient-emphasized image Dr in 

20 which red is emphasized for greater gradients than the gradient image Dra, besides the synthetic 
image Dh of a gray tone indication by a multiplication-synthesis of the aboveground opening 
layer Dp (with ridge emphasized white) and the underground opening layer Dq (with bottom 
emphasized black). 

Then, the gradient-emphasized image Dr is synthesized with the synthetic image Dh. 

25 Fig. 23 is a stereoscopic map of a whole neighborhood of Aokigahara, whereto a 

process according to this mode of embodiment is employed. As seen from Fig. 23, just in the 
south of Tenjinnyama ski area, a glacial hole crater row extends, which has a crater that had 
outflown an Aokigahara lava flow. It is difficult to verify the location in an aerial photograph 
where the deep forest constitutes an obstacle. Further, the glacial hole crater row looks as it may 

30 be in a map by laser data (refer to Fig. 25), while an indication thereof is difficult in a contour 
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map by aerial survey (refer to Fig. 24). 

Contrary thereto, in a stereoscopic picture according to this mode of embodiment, as 
seen from Fig. 6, the glacial hole crater row appears as it is clear defined, as well as a lava flow 
formed with undulations and routes of the climb. 
5 Further, in Fig. 27 and Fig. 28 as magnified views, there can be seen visually defined 

lava flows, road inclinations, and undulations, 

It is noted that the technique in this mode of embodiment is applicable to a terrain of 
Venus as well as a terrain of Mars. In addition, it is applicable for visualization of a ruggedness 
measured by an electron microscope. Further, if applied to a game machine, it allows a 
10 stereoscopic feeling to be given without putting glasses. 

As is above, in this mode of embodiment, based on a DEM (Digital Elevation Model) 
data, three parameters of inclination, aboveground opening, and underground opening are 
determined, and their distribution on a flat plane is stored as a gray scale image. A difference 
image between aboveground and underground is entered in a gray, and the inclination, in a red 
15 channel, for creation of a pseudo-color image, to thereby effect indication to be whitish at a ridge 
or crest part, and blackish in valley or howe, while it gets red as inclined part becomes steep. By 
combination of such indication, an image can be generated with a stereoscopic appearance even 
on a single sheet. Therefore, it allows at a glance a grasp of degrees of concavo-convex heights 
as well as a degree of gradient. 

20 

(Third mode of embodiment) 

Description is now made of a third mode of embodiment of the invention to which the 
visualization processing system VPS1 of Fig. 1 is aided, with reference to Fig. 29 to Fig. 32. 

Fig. 29 is an X-ray fluoroscopic view of a human body, Fig. 30, a graduation-tinted 
25 slice image of the fluoroscopic view of Fig. 29, Fig. 31, a red-toned elevation-depression degree 
distribution image obtained by processing the fluoroscopic view by the visualization processing 
system VPS1, and Fig. 32, a synthetic image of the image of Fig. 30 and the image of Fig. 31 
superposed thereon. 

The graduation-tinted slice image of Fig. 30 is an image, as the X-ray fluoroscopic 
30 view of Fig. 29 having hue of pixels thereof sliced and graduation-tinted in dependence on their 
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brightness, for which, in this mode of embodiment, a field of vectors having positional 
information and brightness of the pixels as their components is stored, in the memory 53, as the 
vector field 70 of the visualization processing system VPS1, and displayed, at the information 
output section 55, as a result of process P7 by the seventh processing file 67 of the visualization 
5 processing program 60. 

Further, the red-toned elevation-depression degree distribution image of Fig. 31 is 
displayed, at the information output section 55, as a result of process P5 by the filth processing 
file 65 of the visualization processing program 60. 

Then, the synthetic image of Fig. 32 is displayed, at the information output section 55, 
10 as a result of process P8 by the eighth processing file 68 of the visualization processing program 
60. 



Industrial Applicability 

According to the invention, there can be provided a visualization processing system, a 
15 visualization processing method, and a visualization processing program, which are adapted to 
visualize a vector field, with local solid attributes thereof inclusive, on a substantially 
two-dimensional field of view, in a manner that allows an intuitive visible perception. 

Further, there can be provided a visualization processing system, a visualization 
processing method, and a visualization processing program, which are adapted to generate a 
20 gradient reddening stereoscopic image that allows at a glance a stereoscopic grasp of terrain's 
heights and inclination degrees 



